Draft version April 26, 2011 

Preprint typeset using I^'T^]X style emulatcapj v. 11/10/09 



COMPONENTS OF THE EXTRAGALACTIC GAMMA RAY BACKGROUND 
Floyd W. Stecker^ and Tonia M. Venters^ 

Draft version April 26, 2011 



o 

(N 
u 

< 
in 

(N 

Oh 
I 

O ' 

(N 
> 
oo 
l> 
\o 
m 



ABSTRACT 

We present new theoretical estimates of the relative contributions of unresolved blazars and star- 
forming galaxies to the extragalactic 7-ray background (EGB) and discuss constraints on the contribu- 
tions from alternative mechanisms such as dark matter annihilation and truly diffuse 7-ray production. 
We find that the Fermi source count data do not rule out a scenario in which the EGB is dominated 
by emission from unresolved blazars, though unresolved star-forming galaxies may also contribute sig- 
nificantly to the background, within order-of-magnitude uncertainties. In addition, we find that the 
spectrum of the unresolved star-forming galaxy contribution cannot explain the EGB spectrum found 
by EGRET at energies between 50 and 200 MeV, whereas the spectrum of unresolved FSRQs, when 
accounting for the energy-dependent effects of source confusion, could be consistent with the combined 
spectrum of the low-energy EGRET EGB measurements and the Fermi-LAT EGB measurements. 

Subject headings: Gamma rays: diffuse background - Gamma rays: galaxies - Galaxies: general — 
Galaxies: active 



1. INTRODUCTION 

Studies of the extragalactic 7-ray background (EGB) 
can provide insight into high energy processes in the uni- 
verse and, as such, has been the subject of much debate, 
particularly concerning the roles of extragalactic astro- 
physical sources and new physics. Recent data from 
the Large Area Telescope (LATj3 on board the Fermi 
Gamma Ray Space Telescope allow for a reassessment 
of the possible astrophysical origins of the EGB, which 
could improve our understanding of 7-ray production in 
these objects and provide more robust constraints on 
the more exotic scenarios. However, in order to deter- 
mine the strength and spectrum of this isotropic back- 
ground one needs to proceed from the raw photon count 
data by determining to the best extent possible the de- 
tector sensitivity, the intrinsic events produced by the 
larger charged particle flux impinging on the detector, 
and the much larger 7-ray foreground within our Galaxy 
resulting from cosmic-ray interactions with photons and 
gas nuclei. Such analyses have been made for both the 
Energetic 7-ray Experiment Telescope (EGRET) aboard 
the Compton Gamma Ra y Observato ry ( Sreckumar et al.l 
ll998tlStrong et al.ll2004aD and Fermi (|Abdo et al.ll2010in . 

Various extragalactic 7-ray production scenar- 
ios have been explored theoretically as candidate 
components that could contribute significantly to 
the observed background. Among those consid- 
ered are the unresolved astronomical sources, such 
as a ctiv e galactic n ucl ei (A G Nl ([Pado vani et al.| 
IQQa IStecker eFall [19931: ISalamon fc Steckefl 
I99I; 'Chiane et al.' '1995^ 'Stecker & Salamon 1993 
Kazanas & Pcrlman 1997: Chiang & Mukhcrieg 
1998; Mukhcrjcc & Chiang 1999; Miickc & Pohl 200a 
Giommi et al.i i2006. : .Narumoto fc Totani 2006; Dermen 



^ Astrophysics Science Division, NASA Goddard Space Flight 
Center, Greenbelt, MD 20771 

^ NASA Postdoctoral Program Fellow, Goddard Space Flight 
Center, Greenbelt, MD 20771 

^ Hereafter we shall refer to the Fermi-LAT instrument simply 
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20071: iPavlidou fc Venteri [20081: Unoue &: Tot anf "2001 
Venters et a l.l 120091: iVentersI 120101: lAbazaiian et all 
20101). star- f qrmiii g galaxies (IPavlidou fc Fields! 120021: 
Fields et all |2010'; 'Makiva et al.' '20111), and star 



burst galaxie s (Thompson et al. 2007; iSteckej 120071 : 
iMakiva et al.l [2011i ). The large majority of associ- 
ated extragalactic sources thus far detected by both 
EGRET ari d Ferrn i are blazars l|Hartmann et al.lll999l : 
lAbdo et all I2010J) . i.e., those AGN for which the 
j et is closely aligned with the observer's line-of-sight 
(jBlandford fc Konigl 119791. including 7-ray loud flat 
spectrum radio quasars (FSRQs) and BL Lacertae-type 
objects. It is expected that since blazars comprise the 
largest class of identified extragalactic 7-ray sources, 
unresolved blazars should contribute significantly to the 
EGB. Additionally, just as our Galaxy produces 7-rays, 
it is expected that 7-rays are produced in other galaxies, 
and as such, unresolved galaxies might also contribute 
to the EGB with the most significant contribution 
originatin g from the population of actively star-f o rming 



galaxies (IStecked [19751 fPavlidoti fc Fieldsl [2OOII [2002 
[Fields et all l2O10i : IMakiva et all 1201 1[ ). Interesting 
truly diffuse mechanisms that could contribute to the 
EGB involve cosmic ray interactions with i ntergalactic 
gas and the cosmic background radiation ( Fazio et all 
[19661: IStecke^ [l97l [DmI [20071: IKeshet et al.l I2003D and 
electromagnetic cascades produced by interactions 
of very high and ultrahigh energy particles with the 
extragalactic backgro und light (Kalash ev et al. . 20091 : 
iBerezinskv eTaD[20lTI : fThlers et al.ll20T6l : lVentersi[2O10f l. 
as well as more exotic scenarios such as dark mat- 
ter annihilation (iSilk fc Srcdnicki 1984; St ecker et al.l 
1985"; 'Rudaz fc Stecker' 198 8: iStecker fc Tvlkal Il989allbl : 
Rudazfc Stecker 1991; U Uio et all 120021) and dec ay 
((Olive'fc Silkl[l985l : iSteckeiilTgM llbarra fc Tran[[2008lf[ . 

In this paper, we estimate the contributions to the 
EGB from unresolved extragalactic 7-ray sources of var- 
ious types and compare them with the EGB obtained 



F or re views on dark m atter annihilation, 
pMB ) and lBertone fc SilM lf2005 ). 



see [Jungman 



et al.l 
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Fig. 1. — Source count distributions for Fermi (black triangles; 
gbdo ct al. 20101) and EGRET (red circles; IReimer fc ThompsonI 
I2001i ) blazars. 



from analysis of Fermi data. In doing so, we also take 
into consideration the effects of both the completeness of 
the Fermi flux limited blazar survey and the important 
effect of source confusion owing to the energy dependent 
angular resolution of the Fermi-LAT detector. We will 
then briefly discuss the implications of possible truly dif- 
fuse emission mechanisms to the EGB. 

2. EGRET AND FERMI RESOLVED SOURCES AND THE 
7-ray LOG N - LOG S RELATION 

In Figure [T] we plot the number of blazars observed 
per square degree versus b lazar flux integrated above 
100 MeV fo r both EGRET (IReimer fc Thomp son 2001) 
and Fermi (|Abdo et al.l I2010J ). In the case of Fermi, 
the source spectra were extrapolate d from a power-law 
fit above a fiducial energy of 1 GeV (|Abdo et al.ll2010k[ ). 
The offset between the resolved source count data ob- 
tained by the two detectors is a result of differences 
in their sensitivit y calibrations ( Thompson et al. 1993,' 
lAbdo et al.ll2010l ). The fiuxes of the Fermi sources are 
reported to be systematically brighter than those re- 
ported by EGRET. 

I 



where Ie{Eq) is the EGB intensity given in units of pho- 
tons cm~^ s~-^ sr~^ GeV~^, L^^max depends on the red- 
shift, source spectrum, and the detector sensitivity, and 
we have differentiated the comoving volume with respect 
to redshift and solid angle, £7. However, as the different 
types of likely contributing 7-ray sources have distinct 
spectral characteristics and redshift and luminosity dis- 
tributions, we must consider each case separately. 

Note that in Equation [T] we have neglected the ef- 
fect of 7-ray absorption due to pair-production interac- 
tions with the extragalactic UV photons. The inclusion 



We note that the 50% completeness in the Fermi sur- 
vey is estimated to be re ached at an effectiv e flux limit 
of - 2 X IQ-^ cm-2 s"^ (|Abdo et al.l 120 lOkD . The flux 
limit varies with galactic latitude and source hardness. 
Nevertheless, the effective flux limit indicated by both 
the tur nover in Fig u re [T] an d the Monte Carlo simula- 
tion in lAbdo etaLl (|2010kD is ~ 2 x 10"** cm'^ s'^ 
By assuming that the Fermi survey is 10 0% complete at 

7 X 10"^ cm~2 s~^ (|Abdo et al.ll2010kD and comparing 
the data for EGRET and the Fermi shown in Figure [1] 
we estimate that the EGRET survey is 50% complete at 
~ 8 X lO"'* cm-2 s"^ 

One may note that the flux level attained for 50% com- 
pleteness for the Fermi extragalactic blazar survey is only 
about four times fainter than that for the EGRET sur- 
vey even though the Fermi-LAT is sensitive to sources 
^ 30 times fainter than that of EGRET. It is important 
to note that there are several factors that determine the 
efficiency of a 7-ray telescope for detecting extragalactic 
sources, among which are: (1) the fiux of the source, (2) 
the spectral index of the source, (3) the intrinsic detec- 
tor background from cosmic-ray induced events, (4) the 
foreground from the Milky Way, and (5) the diffuse ex- 
tragalactic background. The Fermi-LAT was designed to 
reach its optimal effective area for 7-rays with energies 
near and above 1 GeV, whereas EGRET was designed to 
reach its optimal effective area for 7-rays with energies 
near and above 100 MeV. As such, Fermi is more sen- 
sitive to sources with hard spectral indices, particularly 
for the faintest sources observable by Fermi. Thus, it is 
difficult to make a direct comparison between EGRET 
and Fermi. However, we note that the positions of the 
turnovers in the data presented in Figure [1] provide a 
good indication of the observational situation. 

3. DETERMINATION OF THE EGB FROM POINT 
SOURCES 

In general, the total contribution of a given population 
of sources to the EGB is found by convolving the source 
spectrum, Fph(Fo, 2, F^), measured at energy, Fq, of one 
source of 7-ray luminosity, F^, at a given redshift, z, with 
the comoving number density of sources at that redshift 
per luminosity interval, n(F^,z) = cPN/dL^dVcom- We 
then integrate over the comoving volume, V^om, and 7- 
ray luminosity: 



(1) 

I 

of 7-ray absorption would result in a steepening in the 
collective point source spectrum at the hig h-energy end 
(jSalamon fc SteckedlTOQl iVenters erall [2009). though a 
possible contribution from electromagnetic cascade pho- 
tons might m itigate the steepening at higher energies 
(jVentersI l2010f ). However, recent Fermi constraints on 
the 7-ray opacit y imply that the UV background is likely 
to be fairly low ()Abdo et al.ll2010^ . and as such, the ab- 
sorption and the resulting cascades will only have a small 
effect on the EGB. 



Ie{Eo) ^ I I FpY,{Eo,z,L^)n{L^,z)-j-^dL^dz , 
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3.1. The Contribution to the EGB from Unresolved 
Blazars 

In determining the blazar contr ibution to the EGB , 
we foUow the procedure outlined in I Venters et al.l ()2009D . 
We approximate the blazar 7-ray spectrum as a power 
law in energy defined by the photon spectral index, T 
(Fph C!C E~^). The spectral indices of the population of 
blazars form a distribution wi th some spr ead (the spec- 
tral index distribution, SID; iStecker fcSa lamon 1991; 
iVenters fc Pavlidoull2007[ ): hence, the number density of 
blazars is defined as 



n{Lj, z, r) = PjiLy, z)pLiT) = 



dLydVcomd^ ' 



(2) 



where p^{Lry, z) — d^N/dL^dVcom is the blazar 7-ray lu- 
minosity function (GLF) giving the comoving number 
density of blazars per luminosity interval and Pl{^) = 
dN / dV is the normalized blazar SID accounting for spec- 
tral bias (see Section is the 7-ray luminosity 
at the fiducial energy, Ef (taken to be 100 MeV), de- 
fined as E'j times the differential photon luminosity, 
L-ph = d'^N-f/dtdE measured at Ef. is related to 



the integral fiux greater than Ef, F{> Ef), by 

^ AT:D^{r ~ + zfEfF{> Ef), (3) 

where D is the distance measure for the Friedman- 
Robertson- Walker cosmologjO: 

D = -^ r[^A + r!™(l + z')'] dz' . (4) 

A given blazar of 7-ray luminosity, Lj, at a redshift, z, 
with a 7-ray photon spectral index, T, has a measured 
photon fiux of 



FphiEo, z,Ly,r) 



4TTEj[dL{z)]^ 



(5) 

where (iz,(z) = D(l + z) is the luminosity distance. Thus, 
the total contribution to the EGB at a given energy, Eq, 
from unresolved blazars (the collective unresolved blazar 
intensity) is determined by integrating the contribution 
from each individual blazar fainter than the detector sen- 
sitivity. 



lEiEo) 



00 Jo 



Fpi,{Eo,z,Ly, r)p^pL{T)- 



— — -— dL~t dz dr . 
dzdn ^ 



(6) 



where z„ 



5.0 and 



is determined from Equa- 



tion [3] taking F{> Ef) — F^i^, with F^i^ being the 
minimum fiux capable of being resolved by Fermi. We 
have also integrated Equation [T] over the blazar spectral 
index. 

3.2. The Contribution to the ECB from Unresolved 
Star-forming Galaxies 

By applying the same procedure that we use to deter- 
mine the background from unresolved blazars (see Sec- 
tion [33J, we calculate the contribution from unresolved 
star-forming galaxies by determining the 7-ray photon 
flux as a function of energy for one galaxy and then con- 
volving with and integrating over the appropriate cos- 
mological distributions. Is is expected that, as in the 
Milky Way, the 7-ray emission for a star-forming galaxy 
comes mainly from the decay of tt" mesons p roduced 
by co smic-ray interactions with interstellar gas (jSteckeri 
ITOTOt) . The resuhing 7-ray prod uction spectr um has 
been calculat ed by ma ny authors (Stecker' '197i 
1971 [ Cavallo fc Gouldir i971: Stephens & Badhwa! 
Dermcr 1986: Mori '1997; Strong ct al. 2004b 

2010; Kcl ner et al. 2006 : Kama c et al.,.200 6 

,0 



1973 



1981 



, 2001 

Mori 200^. 

For our ca, lc ulatio n we adopt the tt" emissivity given 
by iStecked (|1979| ) renormalized upwards by 25% to 
be co nsistent with the l ocal emissivity measured by 
Fermi lAbdo et al.l (|2009bD . We note that there is also 
emission arising from electron bremsstrahlung, but the 
contri bution is likely to b e small, particularly above 100 
MeV (|Abdo et al.l l2009b[ ). We also neglect the emis- 



5 We take Hq = 70 km s"! Mpc"!, Q„ 

n,. < 1. 



0.3, Ua = 0.7, and 



sion from Compton interactions that may contribute sig- 
nificantly to galaxy sp ectra above 10 GeV, particularly 
for starburst galaxies (iStecked 119771: iHunter et al.l 119971 : 
lAbdo et al.l[2009(] : IStrong et al.ll2010ll . We also note that 
the cosmic ray spectrum should, in fact, vary from galaxy 
to galaxy since it depends on energy dependent leakage 
and each galaxy has a different morphology and mag- 
netic field configuration. This uncertainty can affect the 
predicted slope of the EGB spectrum at high energies, 
but it does not affect the absolute value of the predicted 
background at ~ 200 MeV. 

The 7-ray photon luminosity is related to the 7-ray 
production spectrum per hydrogen atom, q}i{E), by 



Lph(i;) = {qKiE))Nu, 



(7) 



where {qii{E)) is found by averaging q}i{E), the differen- 
tial 7-ray production spectrum per hydrogen atom, over 
the galaxy, and A^h is the number of hydrogen atoms in 
the galaxy in the form of both atomic (HI) and molecular 
(H2) hydrogen. Thus, the 7-ray photon flux for a galaxy 
at redshift z as observed at energy En is 



Fph{Eo,z) 



1 



{qYi[En{l + z)]) Nniz) . (8) 



The rate of production of 7-rays from tt*^ decay is propor- 
tional to the fiux of cosmic rays, which we assume to be 
proportional to the supernova rate. The supernova rate 
is expected to be proportional to the rate of formation of 
higher mass stars, which is, in turn, proportional to the 
overall star formation rate assuming a universal initial 
mass function. 

Assuming then that the rate of production of 7-rays 
from TT*' decay is proportional to the star formation rate 
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(SFR) for a galaxy, we can relate the (gn) of the galaxy 
to that of the Milky Way, (^h^): 



v&(z) 



(9) 



where is the SFR of the galaxy, and we take ((?h^) 
to be some fraction, fq, of the locally measuretfl Qh^- 
We determine fq by integrating the radial profile of the 
flux of cosmic rays weighted by r^. Using the radial 
profiles calculated bv lStecker fc Joiiesl (|1977| ). we obtain 
fq ^ 0.825. Thus, the galaxy spectrum becomes 



The calculation of the amount of gas in a galaxy is 
subject to a considerable degree of uncertainty, especially 
at high redshifts (see Section IT^ . As such, rather than 
focusing on one particular model, we calculate the star- 
forming galaxy contribution for three different models 
arising from different sets of assumptions. In so doing, 
we seek to explore various possibilities and highlight the 
uncertainty. 

3.2.1. Galaxy Contribution Determined from the Schecter 
Function and An Evolving Gas Fraction 

One method for determining the number of hydrogen 
atoms in a galaxy is to assume that the mass of gas in 
the galaxy is some fraction of its stellar mass: 



1 - /gas (2;) mn 



(11) 



0.9 



where fg^z) = A4as/( Mgas + M^) ex ( 1 + z ) 
the gas fraction given in iPapovich et al.l (|2010[ ). and 
we have neglected the possible contribution of helium 
to the gas mass of a galaxy. Substituting Equation 
pT|) into Equation pUj) and approximating vI'(z)/vl/(o) ~ 
Psfr(2:)/psfr(0), we get: 



Fpi,{Eo,z,M^ 



fqq^'^[Eo{l + z)]psFR{z) /gas( 



PSFr(O) 1 - fga,s{z) 

(12) 



where psfr(z) is the cosmic SFR (CSFR) given by 
log(psFR(^)) = -2.06 + 3.391og(l + z) for z < 1 3 and 
PSFRiz) ^ const, for 1.3 < 2: < 4.0 (|Lv et al.ll20'Tll ). 

To get the total contribution to the EGB, we convolve 
with the comoving number density of star-forming galax- 
ies per stellar mass interval as a function of redshift and 
integrate: 



If\Eo)^ I dz^^-^l dM'F^^{Eo,z,M')^z,M'), 
Jo dzdil Jm' . 

(13) 

where $(z,M') = d'^N/dM'dVcom is the Schecter func- 
tion for stellar mass with parameters as determined in 
lElsner et al] (poM ). M' is given by M, = IO^^'A/q, and 

12.0. 



M'. 



we take M^^^ 



8.0 and M^^^ 



^ That is, since is calculated in the literature assuming the 

cosmic ray flux as measured in the solar neighborhood. 

3.2.2. Galaxy Contribution Determined from IR Luminosity 

Functions 

Alternatively, we can determine the 7-ray spectrum 
of a galaxy by assuming that the 7-ray luminosity 
of th e galaxy is prop ortional to so me p ower of its 
SFR, ([Fields et al.l[2()Tnl: IMakiva et al.ll2011l : lAbdo et al.l 
IMoi: iph oc Since Lpi,{E) = {qu{E)) Nu and 

{q-ii{E)) oc 5* (as demonstrated in Section r3. 2. 111 . 



MW 



J{q^^iE))dE^ 



(14) 



where A and a are the best-fit parameters of the above 
power law determined from Fermi observations of star- 
forming galaxies in the Local Gr oup and their SFRs 
(see Section r4.2.2|) . Assuming the iChabrierl (2003) ini- 
tial mass function, the SFR of a galaxy is related to its 
total infrared luminosity, iiR, 



LiR = 1.1 X 10^% 



Mq yr- 



(15) 



ms "eFall 'MO^. The 7-ray photon fiux for a galaxy 
at redshift z is given by 



Fph{Eo,z, Xir) — 



A 



AnD^ \Jq^^dE 



^IR 



q^'^iEoil + z)]. 



(16) 



Then, the total contribution to the EGB is found by 
convolving the galaxy photon flux with an infrared lu- 
minosity function, $(z,Lir) = d^ N / dLmdVcom and in- 
tegrating over infrared luminosity and redshift: 

If\Eo) = J Tz j YQFp^{Eo, z, Lir)$(z, Lir) , 

where we take LiR^niin = lO^^i© and iiR,max — IO^^Lq. 

3.2.3. Galaxy Contribution Determined from the Cosmic 
Star-formation Rate and the Star-formation 
Efficiency 



Another alternative is to relate the cosmic density of 
hydrogen in star-forming galaxies to the cosmic star for- 
mation rate. Given that stars are formed in giant molec- 
ular clouds (GMCs), it is reasonable to assume 

PSFR ^ C(H2)PH2 , (18) 

where psfk is the cosmic star formation rate den- 
sity (see Section I3.2.ip . is the cosmic molecular 
hydrogen density in star- forming galaxies, and C(H2) 
is the star formation "efficiency" (SE E) of molecu- 
lar hydrogen ( Bigicl et al. 2008; Gne din et all 120091: 
IBauermeister et al.ll2010l) . iLerov et al.l (j2008[) measure 
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the SFE to be - (5.25 ± 2.5) x 10"i"yr-i and to be 
roughly constant over a wide range of condition^ We 
can relate the density of atomic hydrogen to the density 
of molecular hydrogen through the average mass ratio 
of atomic and molecular hydrogen (7^ = (Mhi/Mhj)) in 
star-forming galaxies, pm ^ TZpii^. The average mass 
ratio of atomic and molecular hydrogen can be found by 



J 



integrating radial profiles of the gas surface den sities of 
star-forming galaxies found in lLerov et al.l (|2008D . result- 
ing in TZ ~ 0.9. Note that in so doing, we only integrate 
the profiles out to the optical radius since recent sur- 
veys indicate that star form ation is ex tremely inefficient 
beyond this radius (Bigi el et~ani2010f fl. Thus, with ap- 
propriate modifications to Eguation llOl we find the 7-ray 
flux from a particular redshift 



PsFr(^) dVco, 

Psfr(O) dz 



-dz . 



(19) 



r 



Differentiating with respect to solid angle £7 and inte- 
grating over redshift results in an equation for the total 
contribution to the EGB: 



47rmH^(H2)/3sFR(0) 



-1 

ZJ2 



,9lf^[i?0(l + z)]piFRW 



dfldz 



■dz . 



(20) 



r 



4. OBSERVATIONAL INPUTS AND CONSIDERATIONS 

4.1. ^-ray Blazars 

Of the nearly 1500 resolved point sources observed 
by Fermi in the first year, 573 a re associated wit h 
blazars (First Fermi catalog flFGL): lAbdo et al.|[20lM ). 
Thus, blazars comprise the largest class of astrophysical 
objects associated with 7-ray sources. Naturally, unre- 
solved blazars have long been suspected of providing, 
at least, a substantial contribution to the EGB, though 
the exact amount remains in debate and depends on 
various assumptions as to c onstruc ting GLF s and red- 
shift d istributions (Padova ni et al.l |l993: St ecker et al.l 

IChiang et al.l 11995 



19931 iSalamon fc Steckeri 119941: IChiang et"alT 

Stecker fc Salamon 19961 Kazanas fc PerlmanI 
^ — — 1^ — — ^^^^^^^"^^^^^^^^ 



Chiang & Mukheriee 
Mukhcricc & Chiang 



1999: 



'Sreckum ar et aLl 
Miicke fc Pohj 



1997 



1998 



2000; 



Giommi et al.l 12 006: Narumoto fc Totani 2006; Dormer 
20071: iKneiske fc Mannheimi i2008: Pavlidou fc Venters 
200a llnoue fc Totanil 120091: \ ^nters et al.l 120091: 
Abdo et all l2010kl : iVentersI I2010D T A detailed dis- 



cussion of all of the assumptions that go into these 
calculations is beyond the scope o f this paper (though, 
for a detailed discus sion of the IChiang fc Mukheried 
119981 calculation, see iStecker fc SalamonI l2001| ). and it 
is likely that many of these models will be updated 

More precisely, 5(H2) ~ e/Tfl-, where e is the percentage of 
gas involved in forming stars and oc p~^/'^ is the local free-fall 
timescale of the gas. However, using this relation requires knowl- 
edge of the local density of the gas and a better underst andi ng of 
the formation of GMCs than presently exists (see Section 14.211 . We 
also note that the measurements obtained by[Lerov et al. (200^) 
were taken from a sample of low-redshift galaxies. The SFE could 
actually evolve with redshift (Bauermeister et al. 2010). 

* We should note that even though star formation is extremely 
inefficient beyond the isophotal radius, there is still gas beyond 
this radius. However based on the radial profiles determined in 
[Stecker & Jones (1977), we do not expect cosmic rays to propagate 
much beyond the isophotal radius; hence, we expect gamma-ray 
production beyond this radius to be low. 



in light of Fermi data. Howev er, we note that as 
discussed in |Abdg_et_alJ |2010kD . the source counts 



3% , 

predicted by Dermefl |2067|) and iMiicke fc Pohl pOOl 
fall short of the Fermi observations of resolved sources 
above 5 x 10~^ phcm"^ s~^. In calculating the blazar 
contribution to the EGB, we assume functional forms 
for the blazar GLF and SID and fit them to IFGL data, 
accounting for errors in measurement of blazar spectral 
indices and the spectral bias inherent in a flux-limited 
catalog. 

4.1.1. Source Counts for Faint Unresolved Blazars: Theory 
Meets Observations 

Determining the GLF from observations relies on the 
ability to associate 7-ray blazars with lower energy coun- 
terparts f or which redshifts can be measured (for discus- 
sion, see iVenters et al.l 120091 : IVentersI 120101 ). However, 
making the necessary association can be complicated by 
the angular resolution of the Fermi-LAT, which is limited 
by electron scattering in the LAT detector and is much 
poorer than that of more traditional telescopes (see Fig- 
ure [3]). The resulting wide point-spread function results 
in significant source confusion at energies below 1 GeV 
even for fluxes well above the Fermz-LAT sensitivity. In 
principle, one could construct source counts from fluxes 
integrated above an energy for which source confusion 
is less of a hinderance, but doing so would limit the al- 
ready suppressed blazar number statistics. Thus, rather 
than construct a luminosity functio n solely from 7-ray 
blazars with redshifts, we employ the IStecker fc SalamonI 
(|1996f PI approach of determining the luminosity function 
from wavebands with larger samples and smaller posi- 
tional er ror circles. 

As in IStecker fc SalamonI (|1996f ). we take the func- 
tional form o f the FSRQ luminosity fu nction from radio 
observations (|Dunlop fc Peacockll990( ) , but corrected for 

^ See also lNarumoto fc TotaiTil 1(2009 ). 
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the present cosmological parameters. The 7-ray luminos- 
ity of a blazar is then determined from its radio luminos- 
ity. The average correlation between the radio and 7-ray 
luminosities of blazars is determined by fitting the bright 
end of modeled source count^ to that of the observed 7- 
ray source counts (Xreducod ~ 0.4). In so doing, we find 
that 7-ray luminosity integrated from 100 MeV to 100 
GeV is ~ lO'^'^ times vL^ in radio, in agreement with 
the results obtained using recent Fer mi observational 
LGiroletti et al."2010'; 'Abdo et al.ll2010il : iGhirlanda et all 
[2010; Mahony ct al. 2010). We should also note that 
we include sourc e s out to 2; ~ 5, but as demonstrated 
in iVenters et al.l (|2009f) . the emission is dominated by 
sources with z < 2 (consistent with expectations from 
the redshift distribution predicted by the GLF). As such, 
we do not expect this choice to significantly impact the 
results. 

The resulting modeled source count distribution for 
FSRQs (sohd) is presented in Figure [5] along with the dis- 
tributions of the brighter FSRQs resolved by Fermi (light 
data points) and all blazars (dark data points). We also 
show the forms of the unresolved source count distribu- 
tions obtained from a Monte Carlo modeled Jer 7 7t^-LAT 
sensitivity calculation performed bv lAbdo et al.l (j2010kD 
as dashed lines. All of the models separate from the data 
at fainter fluxes as the source counts fall off very rapidly 
due to the survey incompleteness and source confusion. 
Thus, the determination of the true faint-end shape di- 
rectly from the source counts can be severely hindered. 

In an effort to mitigate the effect of survey complete- 
ness, lAbdo et all (|2010k[ ) modeled the Fermi-hkT effi- 
ciency from a Monte Carlo simulation and then divided 
the differential source counts by this efficiency. Thus, the 
calculated source counts for faint sources strongly depend 
on such modeled efficiencies. Indeed at fluxes near the 
Fermi-LAT sensitivity limit, the efficiency is extremely 
small and model dependent; hence, in flux bins with low 
number statistics, the source counts are multiplied by a 
very large and uncertain number. The result of such a 
procedure is that even though one source is not statis- 
tically different from two sources, whether one source is 
seen or two could result in different modeled counts for 
faint sources. It has also been argued that unassociated 
7-ray sources are likely to be domina ted by known classes 
of 7-ray sources, especially blazars d Mirabal et al.|[20To() 
and will have a contribution to the EGB even though 
they will not have been included in the source count dis- 
tributions. It is also important to note that in determin- 
ing the blazar contribution to the EGB from measured 

^'^ In doing so, we also include the blazar SID (see Section [4. 1.31 1. 
One might be concerned that the effect of the new radio-7 
correlation would be to increase the blazar background with re- 
spec t to the Stecker & Salamon (1996) model. However, we note 
that [Stecker & Salamon, (1996.) distinguished between "quiescent" 
and "flaring" blazars and used separate radio-7 correlations (and 
spectral properties) for each subpopulation. In this paper, we 
make no such dis tinction, so a com parison between our results 
and those of the Steck er &: SalamonI (11996) model is not straight- 
forward. Using the Stecker & Salamon (1996) radio-7 correlation 
for quiescent blazars from Stecker & Salamon ( 1996) would yield 
a smaller blazar background, but it would also under-predict the 
bright end of the 7-ray source counts. To compensate, on e woul d 
have to add a flaring component (as per lStecker fc Salamon! ( 119961) ) 
to fit the data, which would also contribute to the background. In 
effect, such a procedure is equivalent to our method of fitting the 
radio-7 correlation to the data. 



source count distributions, source confusion could not be 
taken into account since its effect depends on the source 
density (see Section l4.1.2[) . which is exactly the unknown 
quantity that the observer seeks to determine. As such, 
it is likely that the blazar contribution to the EGB will 
be underestimated in analyses based solely on the mea- 
sured source count distributions. Since in this paper we 
use a theoretically determined source count distribution, 
our model gives a source count density from which one 
can determine the effect of source confusion (see Section 
14.1.21) . Sucli differences bet ween our analysis and that 
given in P\_bdo et al.l (|2010k[ ) result in different calcula- 
tions of the unresol ved blazar contribution to the EGB: 
I Abdo et"al] (|2010k[ ) conclude that blazars can only ac- 
count for less than 25% of the EGlj^. while our analysis 
indicates that blazars could possibly account for the bulk 
of the EGB. 

The -Fioo fluxes included in the source count distri- 
butions presented in Figure [1] are not actually measured 
i^ioo fluxes since the IFGL catalog does not include i^ioo 
fluxes. In order t o determine t he J"ioo fluxes, we make use 
of Equation 1 of lAbdo et al.l ()2010kf) to extrapolate Fioo 
fluxes from mea sured differe n tial flu xes determined at the 
pivot energies. lAbdo etall ()2010J ) deflne the pivot en- 
ergy as that energy for which the differential flux is min- 
imal. We note that for IFGL blazars, the average of the 
pivot energies is ~ 1 GeV. Thus, the source count distri- 
butions might be more representative of source counts for 
sources brighter than the Fermi sensitivity above 1 GeV 
rather than 100 MeV. In effect (and also as the result of 
spectral bias; see Section ll.l.3p . the source counts could 
underestimate the number of blazars with i^ioo fluxes 
above the Fermi sensitivity, mostly impacting the faint 
end of the source count distributions. As such, analy- 
ses on unresolved blazars based solely on source counts 
likely underestimate the blazar contribution to the EGB. 
This effect could provide an explanation for the fact that 
Fermi observes roughly as many BL Lacs as FSRQs even 
though BL Lacs are intrinsically fainter than FSRQs. 
Since the 7-ray spectra of BL Lacs are harder than those 
of FSRQs they are more easily observed by Fermi than 
FSRQs. 

4.1.2. The Fermi-LAT and EGRET Angular Resolutions 
and Source Confusion 

As mentioned previously, the large angular resolutions 
of pair-production 7-ray detectors such as EGRET and 
the Ferm«-LAT0 result in significant source confusion, 
particularly for faint sources and for energies below ^ 1 
GeV. Since the angular resolution for the Fermi-LKT is 
similar to that of EGRET at - 100-200 MeV, the capa- 
bility of Fermi to resolve faint sources is similar to that of 
EGRET at these energies. Hence, if the EGB does indeed 
consist of unresolved sources, then the EGB measure- 
ments of the detectors should be similar at ^ 100 — 200 
MeV, whereas at ~ 1 GeV, the improved angular reso- 
lution of Fermi with respect to that of EGRET should 
result in a lower measurement of the EGB by Fermi than 
that of EGRET owing to the enhanced ability of Fermi 

12 In Section 8 of Abdo et al.' f2010k) a value of 40% of the EGB 
is obtained if one extrapolates to zero flux. 

,ht tp : / /www-glast . slac . Stanford. edu/software/IS/glast_| 
|lat_perf ormance .htm . 
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Fig. 2. — Bright end of the source count distributions for all 
blazars (including BL Lacs; black data points) and FSRQs (light 
green data points). Our model fit to the data is shown by the solid 
(purp le) line. The dashed lines are the faint-end slopes determined 
in A bdo et al.l 1120101a ') by including a modeled Monte Carlo Fermi- 
LAT efficiency. 



to resolve point sources (jStecker fc SalamonI [19991) . 

The probability, V, of finding a nearest neighboring 
source with S > Sum within the minimum angular sepa- 
ration, ^min, for a source density, N , is given by 

r{< OnuniE)) = 1 - exp(-7riVCi„(i?)). (21) 

For our source confusion criterion, we take the acceptable 
probability limit, "Pmin ~ 0.1, and 6'min(£^) 6(,j%{E) 
approximately given by 



eQ^%{E) = 5.12° X 



E \ 



100 MeV 



(22) 



(|Atwood et all [20091 ) ( see Figure [3|) . Then, the source 
density criterion (SDC) is found by inverting Equation 



NsDC 



ln(l - 



(23) 



The limiting source flux, S'sdCj is then determined from 
the modeled source counts. 

If one were to think of the Fermi-hAT as an ordinary 
telescope with = 6*67% where ^67% is the half angle 
for a beam the contains 67% of the photons. The source 
density criterion would correspond to ~ 1/10 sources per 
beam. For sources with Fioo ^ 1 x 10~^ph cm~^ s~^, 
the probability for finding another 7-ray source of sim- 
ilar or greater flux within the error circle is quite high; 
hence, many sources that, in principle, should be resolv- 
able are, in fact, unresolved. At fluxes close to the sensi- 
tivity limit, this probability is so high that faint sources 
are indistinguishable and will contribute to the measured 
EGB of the detector in question (for the Fernii-LAT res- 
olution at 100 MeV, the source criterion corresponds 
to A^sDC 1-3 X 10^"^ sources/deg^; for reference, at 
~ 2 X 10~^ phcm~^ s~^, our model predicts a much 
larger source density of ~ 0.3 sources/ deg^). Thus, the 
effect of source confusion compounded with the separate 
effect of detector sensitivity is to flatten the faint end of 
an observationally derived source count distribution. 




Energy [GeV] 



Fig. 3. — The angular resolu tion as a function of energy for 
the Fermi-LAT ( solid black lineilAtwood et al.ll2009D and EGRET 
(dashed red line; [Thompson et aLHToQ^ 



We should emphasize that this definition o f source con- 
fusion i s not the same as that employed bv lAbdo et al.l 
()2010kl) . Source confusion as discussed in this section 
refers to the probability that for a given source density, 
a source has a nearest neighbor within the angular reso- 
lution of the det ector with a flux gre ater than or equal to 
the flux limit. In lAbdo et al.l (|2010k[ ). the term is applied 
to a detected source associated with a given real source 
for which the measured flux of the detected source is 
greater than the flux of the real source (plus three stan- 
dard deviations) by a given amount (the analysis was 
performed on simulated data to determine the impact 
on the actual data). In effect, the former definition iden- 
tifies the limit at which the source density is sufficiently 
large so that individual sources cannot be resolved, and 
only fluctuations are observed. The latter definition ap- 
plies to the probability that a given detected source is in 
fact the superposit ion of several sources. However, the 
I Abdo et"an (|2010kl ) criterion is based on the assumption 
that a source can be resolved. As such, the source would 
have to be significantly brighter than the background, 
including the sources within the error circle of the de- 
tector. In the case of many sources with similar fiuxes, 
none of the sources would be resolved. In the case of one 
bright source and several fainter sources, the flux of the 
detected source would be dominated by th e flux of the 
brighte st source and likely would escape the lAbdo et"al] 
(l2010k[) . Thus, under this criterion, it is not surprising 
that they conclude that there is very little source confu- 
sion. A common treatment of source confusion in mea- 
sured source counts is a fluctuation analysis, but as yet, 
such an analysis has not be performed on Fermi source 
counts. 



4.1.3. Blazar Spectra 

The spectral indices of the population of blazars 
form a d istribution wit h a given finite sp read 
(jStecker fc Salam on 1996; Venters fc Pavhdoul[2007l) . re- 
sulting in curvature in the spectrum of the collec- 
tive intensity of unresolved blazars due to the increas- 
ing relative importance at high energies of blazars 
with harder sp ectral indices (jStecker fc Salamonl 119961 : 
iPavlidou fc Venters ,20081) . Thus, the determination of 
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the blazar SID is crucial in determining the correct spec- 
trum of the unresolved blazar contribution. 

In determining the blazar SID from survey data, one 
must carefully account for uncertainties in measurement 
of the spectral ind ices. In so doing ^ we follow the likeli- 
hood method of Venters fc Pavlidoul p007), fitting Fermi 
IFGL FSRQs to a gaussian SID. We determined the 
maximum-likelihood gaussian SID parameters (mean, 
To, and spread, ctq) to be Fq = 2.45 and ctq — 0.15. 

We must also account for the effect of spectral bias 
of the IFGL catalog. Even in a flux-limited catalog, 
low-luminosity, high-redshift blazars that would be most 
likely to appear are those with spectral indices that are 
harder than most of the population. The IFGL catalog 
is not actually flux limited (see Section [2]), rather it con- 
sists of sources above a given threshold in test statisticF^. 
which depends on source s pectra and the background . 
However, as demonstrated in lVenters fc Pavlidoul ()2010D . 
applying the likelihood analysis to the sample of FSRQs 
with fluxes > 7 x 10~^ ph cm~^ s~ ^ and galactic latitudes 
> 10° (as per lAbdo et al.l[2010k[ ) does not appreciably 
change the SID parameters (Fq = 2.45, ctq = 0.16^3. 
Thus, even though the IFGL catalog is not flux limited, 
we can assume that the sample of IFGL FSRQs is ap- 
pr oximately flux limited. We can then follow the method 
of iVenters et al.l ()2009[ ) in correcting for the sample bias 
inherent in a flux-limited catalog. In doing so, we apply 
a correction f actor, M{a), to the SID (for derivation, see 
IVenters et alll2009l) : 



PL (a) = 



M{a) 



(24) 



where f>{a) is the SID corrected for measurement uncer- 
tainty in the spectral indices, and 



M{a) oc 



1 



7 J z=a 



azp^{a, z, r^) — - — 



(z). 
(25) 



4.1.4. Summary of Differences with the lStecker & Salamori 
(1996) Blazar Model 

To summarize, our approach in calculating the 
blazar con tribution to th e EGB is similar to that of 
I Stecker fc~S alamonl (|1996D with some notable differences: 

1. The model has been updated to make use of the 
current cosmological parameters. The change in 
cosmology has little impact on the results. 

2. The model considers only FSRQs, and flaring 
blazars are not considered separately from quies- 
cent blazars. (For possible impact, see item 3.) 

3. The 7-ray-radio relation has been updated to be 
consistent with multi-wavelength observations of 
FSRQs conducted by Fermi and radio telescopes 
(L^ ^ 10 ^-^ X vLy as opposed to 10 ^-^, see Sec- 
tion UTTH])- As noted in Section [4.1.11 the smaller 
7-ray-radio relation would result in a model that 
cannot by itself fit the data and would require a 

For IFGL, this threshold is 25, corresponding to a statistical 
significance of ~ 4(7. 

15 See also lAbdo eTaLl H2010kl ). 



fiaring component. In effect, this is, on average, 
equivalent to our choice of a single population of 
blazars with a given 7-ray-radio relation. We do 
not expect this choice to have a major impact on 
our results. However, as we will discuss in Sec- 
tion |4?T|5l the blazar duty cycle is a remaining un- 
certainty, and it could impact predictions for the 
number of blazars that are observable by Fermi. 

4. The model has been updated to account for effect 
of source confusion in the in blazar contribution 
to the EGB (see Sections |4. 1 . 2! and 15]) . As wiU be 
shown in Section |5l this has the effect of increasing 
the blazar background at lower energies. 

5. The SID o f FSRQs has been updated following the 
analysis of IVenters fc Pavlidoul (|2007D (Fq = 2.45, 
ctq = 0.15) and cor recting for spectral bias as in 
IVenters et all (|2009f) (see Section EO]). Thus, the 
collective spe ctrum of blazars is n o t as h ard as that 
presented in iStecker fc SalamonI ()1996l ) and does 
not exhibit as much curvature. 

4.1.5. Remaining Questions 

There remain a few open questions, the answers to 
which will impact the determination of the blazar con- 
tribution to the EGB. The blazar duty cycle, which dic- 
tates the amount of time a blazar spends in the quiescent 
state versus the fiaring state, remains uncertain, as do 
questions of the amount the fiux increases during fiaring 
and whether the spectral index changes during fiaring. 
Analyses of EGRET blazar spectral indices found no ev- 
iden ce of systematic changes in spectral index with fiar - 
ing (jNandikotkur et al.ll2007l : IVenters fc Pavhdoull20Q7| ). 
and Fermi observations of individual blazars have thus 
far revealed no systematic changes in spectral inde x with 
time o r fiux (Abdo et al. 2009d. 2010d: Acker mann et al.l 
HoToB). As such, we are justified in assuming that the 
blazar spectral index remains constant, on average, with 
time. However, we acknowledge that the uncertainty of 
blazar variability parameters could have an impact on 
the counts of faint blazars and the 7-ray-radio correla- 
tion. This uncertainty will decrease as more data from 
Fermi become available. 

Another open question is that of the nature of blazar 
spectra over the entire Fermi energy range. We treat 
blazar spectra as unbroken power laws over this range 
and in many observed blazars, this does appear to be 
a reasonable approximation. However, in at least a few 
cases, Fermi has found evidence that blazar spectra can 
break (|Abdo et al.ll2009al . l20Toil ll). Whether such obser- 
vations are representative of the entire blazar population 
is presently unclear, as is the nature of the breaks. In any 
case, spectral breaks are likely to impact the collective 
unresolved blazar spectrum mostly at the high end of 
the Fermi energy range. It is also possible that the spec- 
tra of harder blazarq3 will compensate for that of so fter 
blazars at higher energies (IVenters fc PavhdoullMol ). 

4.2. Star-forming Galaxies 

1^ We do not include a possible contribution from BL Lacs for 
lack of a comparable radio data set. 
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As discussed in Section [X2l the Milky Way is a source 
of substantial 7-ray emission arising primarily the decay 
of 7r° meson produced in inelastic collisions of cosmic 
rays with interstellar gas. Thus, it is expected that other 
star-forming galaxies emit 7-rays through the same inter- 
actions and that unresolved star-forming galaxies could 
provide a substantial contribution to the EGB. Thus far, 
the Fermi-LAT Collaboration has reported detections of 
two nearby irregu lar galaxies (the SMC and the LMC; 
lAbdo et all l2010clllJ) two starburst galaxies (M82 and 
NGC253: IAbdo et al.l[2010b'). a nd M31, a galaxy similar 
to our own lAbd o et al.ll2010gD . As such, whatever the 
contribution to the EGB from unresolved star-forming 
galaxies, it will not have changed substantially in the 
Fermi data with respect to EGRET data as Fermi has 
resolved only a handful of star-forming galaxies. 

The determination of the star-forming galaxy contri- 
bution relies on knowledge of the star formation rate of 
galaxies and their gas content, both of which are sub- 
ject to substantial observational and theoretical uncer- 
tainties. At relatively low redshifts (z < 1.5), nebular 
and forbidden emission lines (e.g.. Ha, O II, and O III) 
can be used to trace star formation in galaxie^. At 
higher redshifts (z ^ 1 — 5), the redshifted UV contin- 
uum is used. However, both observational techniques are 
subject to uncertainties in dust extinction and the stel- 
lar initial mass function. Alternatively, in noting that 
UV radiation from young stars is absorbed by interstel- 
lar dust and reradiated in the infrared, measurements of 
the far-infrared (FIR) continuum can be used to trace 
star formation. Nevertheless, early-type galaxies can 
exhibit substantial FIR emission possibly due to dust 
heating by older stars or AGNs. Furthermore, infrared 
measurements are hindered by emission from our own 
Galaxy. Given these factors, the large degree of scatter 
present in the measurements (about a factor of a few; 
iLe Borgne eran[2009t iLv et alJiMll ) of the cosmic SFR 
density is not surprising. 

The gas content of galaxies is even more uncertain, 
particularly at high redshift. The amount of H2 in a 
galaxy is determined from measurements of CO emis- 
sion, while HI is determined from measurements of the 
21-cm line. However, the CO-to-H2 conversion varies de- 
pending on the metallicity and radiation field of a given 
region and the opacity of the molecular clouds containing 
CO. The 21-cm surveys, on the other hand, extend only 
out to 2: ~ 0.05. At higher redshifts, measurements of 
the HI density of the universe rely on damped Lyman-a 
absorbers observed in the Lyman-a forest of quasar spec- 
tra, but the nature of these systems is still the subject of 
much debate (see e.o.. lKulkarni et aD2010l : IPeroux et al.l 
l2010f ). Furthermore, the connection bet ween the total 
gas a nd star formation rate is complex (jPutman et al.l 
l2009f ). While it is fairly well established that stars form 
in GMCs and hence that star formation traces H2, the 
amount of HI varies between galaxies and does not ap- 
pear to be correlated w ith star formation (Bigie Fet al.l 
120081: iLerov et al.]|2008[ ). From both the observational 
and the theoretical points of view, the transition from 
HI to H2 and the formation of GMCs remain uncertain 
(|Lerov et al.ll200l . 

^"^ For review of observational techniques of measuring star for- 
mation, sec Kcmiicutt (1998.). 



Given the uncertainty surrounding key elements of the 
determination of the star- forming galaxy contribution to 
the EGB, our approach does not focus on a particular 
model. Instead, we employ several families of models 
that rely on separate sets of assumptions each with ad- 
vantages and caveats. Using this approach we seek to 
explore various possibilities for the star-forming galaxy 
contribution and highlight the uncertainty in such a cal- 
culation. The strategies we employ are summarized as 
follows: 

1. Relate the galaxy gas mass to its stellar mass as- 
suming a gas fraction that evolves with redshift. 

2. Relate the galaxy 7-ray luminosity to its SFR, 
which, in turn, is related to an observable for which 
there is a redshift distribution {e.g., IR luminosity). 

3. Relate the cosmic density of gas in star- forming 
galaxies to the star formation rate density. 

4.2.1. The Schechter Function Model 

In this approach, we relate the galaxy gas mass to its 
stellar mass assuming an evolving gas fraction (for de- 
tails of the model, see Section 13.2. ip . We employ the 
Schechter parame ters of th e stellar mass functions as 
determined by Eis ner et al. | (2008) and th e evolv ing gas 
fraction as determined bv IPapovich eiTall (1201(1 . bmce 
extensive sp e ctrosc opic surveys are, as yet, unavailable, 
I Eisner et al.l ()2008l ) make use of combined data from the 
multi-band photometry of the GOODS-MUSIC catalog 
and the Spitzer Space Telescope. In so doing, they in- 
fer the stellar masses of galaxies from photometric data 
by fitting the mass-to-light ratios of galaxies to stellar 
population templates. Such a procedure is subject to 
a considerable degree of uncertainty, particularly arising 
from that of dus t extinctioiF^ and the usage of photo- 
metric redshifts. lElsner et al.l ()2008[ ) estimate the mean 
uncertainty in their stellar mass estimates to be about 
a factor of two, though they did not estimate the possi- 
ble uncertainty resulting from the usage of photometric 
r edshifts. 

IPapovich et al.l ()2010f ) study the relationship between 
the star formation rate and stellar mass of high redshift 
galaxies selected at a constant comoving number density 
and derive a gas fraction that evolvej^ as (1 -I- z)°-^. 
They estimate that the uncertainty on the gas mass is ^ 
0.11 dex. 

The advantage of this kind of model is that it is based 
on observations that will be continuously refined with 
time. However, we note that the relationship between 
the stellar mass of a galaxy and its total gas content is 
unclear, though perhaps with better observations and a 
better theoretical understanding, the relationship will be 
better determined. 

4.2.2. The IR Luminosity Function Models 

In this approach, we assume that the 7-ray luminos- 
ity of a star-forming galaxy can be related to its star 

In order to mitigate th e effect of the uncertainty in dust ex- 
tinction, |ElsneF^F|al] II2OO8I ) make use of Ks-band M/L ratios since 
dust absorption is small for longer wavelengths. 

Note that we extrapolate the functional form of their gas 
fraction to low redshifts. 
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Fig. 4. — The 7-ray luminosities of Local Group galaxiesplotted 
versus their star formation rates (data points; see Table [TJ. Also 
plotted is the power-law fit relating the 7-ray luminosity to the 
SFR (solid red Une): L-y oc 



TABLE 1 

Observables for Local Group Galaxies 



an IR luminosity function, and integrate with respect to 
IR luminosity and redshift. 

While the use of IR luminosity functions taken from 
observations is possible, it is difficult to deconvolve the 
contributions from obscured AGNs and mergers. As 
such, we use the semi-empirical IR luminosity functions 
determ ined from the halo-o ccupation-based methodol- 
ogy of iHopkins et al.l (|2010[ ) for both star-forming and 
starburst galaxies (in which enhanced star formation 
due to major m ergers is taken into account While 
the IHopkins et a l. (2010) IR luminosity functions match 
available observations fairly well, we note that there is 
considerable debate over the roles of AGNs and mergers 
in driving star formation, the evolution of galaxies, and 
the determination of the IR luminosities of massive sys- 
tems. Furthermore, the 7-ray-SFR correlation is based 
on rather uncertain estimates of the SFRs and 7-ray lu- 
minosities of one normal galaxy (our own), two irregular 
galaxies, and two starburst galaxies, all of which could, 
in principle, exhibit different star formation properties. 
As more data become available from Fermi, the 7-ray- 
SFR correlation will be further tested. If it proves robust, 
then studies of the normal galaxy contribution the EGB 
could have implications for large-scale-structure forma- 
tion and the evolution of galaxies with cosmic time. 





L~i (ph s -"-) 


*(M0yr-i) 


SMC 


(1.7 ±0.3) X 10''" 


0.01 ± (4.8 X 10-*) 


LMC 


(7.8 ± 0.6) X lO^'O 


0.1 ± (4.6 X 10-=*) 


M31 


(6.6 ± 1.4) X 10"'! 


0.4 ±0.18 


Milky Way 


(1.2 ±0.2) X 10"'^ 


0.93 ±0.34 


NGC253 


(1.1 ± 0.7) X 10*3 


2.7 ±0.5 


M82 


(2.5 ± 0.9) X 10''^ 


5.4 ± (5.9 X 10"*) 



^_3MC: Lj from lAbdo et al.l Il2010d 1. 
l[2OT0l 'l corrected for IMF. LMC: L, 



SF R from iLawton et al.l 
from lAbdo et al.l l|2010hl 1. 



SFR from ILawton et al.l l|2010) corre cted for IMF. M31: L. 



from 



Abdo et al.l Il2qi0e l). SFR from IT abatabaei BerkhuTisenI 



H20ld ) and IWilliamsl "12003') corrected for IMF. MW: calcu- 



lated for (qh) ~ 1.4 X lO-^^s-^H-i and A/h - 7 X lO^Afn 
i Boissier &: Prantzosjig^ gg). SFR taken from Robitaillc & Whitney 
I 20101 ) corrected for IMF. NGC253: from Abdo et al. (2010b). 
SFR fro m IR measu rements given by Sanders et al. (2003). M82: 
L-Y from'Abdo et al.J ('^OlOb). SFR from IR measurements given by 
[Sanders et al. (2003). Note that the errors in the SF Rs r eflect only 
statistical uncertainties and assume that Equation 1151 is exactly 
correct. 

formation rate as a power law. In order to determine 
this relationship, we fit 7-ray luminosities of Local Group 
galaxies calculated from Fermi measurement^ to their 
star formation rates either taken from the literature 
or calculated from IR measurements and converted to 
SFR via Equation [ini In all cases, the SFR is calcu- 
lated assuming (or corrected to) the IMF of iChabrieil 
(|2003D . We find the best-fit power law to be given by 
Lj [lO^^hs-i] ~ 24.0 X ^-i-^. This fit is consistent with 
that obtained by th e Fermi Collaboration, L-y oc \[/i-4±o.3 
(|Abdo et al.ll2010^ . Observables for Local Group galax- 
ies used in this analysis are given in Table [1] and plotted 
with the fit in Figure 21 In a manner similar to that 
of blazars, we can relate the 7-ray luminosity of a star- 
forming galaxy to its total IR luminosity, convolve with 

For the Milky Way, we calculate the 7-ray luminosit y fro m gjj 
averaged over the whole galaxy (as discussed in Section 13.21 1 and 
the Afjj taken from IBoissier fc Prantzo3 (119991) . 



4.2.3. The Strong Coupling ^ -ray-Star Formation Rate 

Model 

In this approach, we seek to relate the cosmic den- 
sity of hydrogen in star-forming galaxies to the cosmic 
star formation rate. Observations of nearby galaxies 
have indicated that localized star formation traces the 
density of H? but there is no direct correlation with HI 
(|Bigiel et al.f2008l : iLerov et al.ll2008| ). While in principle, 
there should be some relationship between the amount of 
HI, the amount of H2, and the SFR in a galaxy, other 
factors such as density fluctuations and turbulence make 
such a relationship complex. Therefore, while acknowl- 
edging that the relationship between HI and SFR is un- 
clear, we simply assume that the amount of HI in star- 
forming galaxies is, on average, comparable to that of 
H2 within the optical radius of the galaxy (see Section 
13.2.31) . Hence, in this model, we take the 7-ray lumi- 
nosity to be roughly proportional to the square of the 
SFR. We stress that the assumption that the star for- 
mation rate is proportional to the available gas density 
only applies to galaxies that are actively forming stars. 
This assumption does not apply to galaxies at very high 
redshifts which may contain substantial amounts of gas 
but have yet to begin forming stars. However, we only 
include star-forming galaxies out to z ~ 4, so the higher 
redshift galaxies will not impact on our results. Given 
that the best-fit power-law index determined for Local 
Group galaxies in Section [4.2.2! is ~ 1.2 and given the 
proximity of the resulting unresolved spectrum to the 
Fermi measurements of the EGB (see Section [5]), we con- 
sider this model to be reflective of an upper limit to the 
star- forming galaxy contribution to the EGB. 

■^^ Note that we do not consider any contribution from the so- 
called calorimetry effect for starburst galaxies as such an effect is 
likely to be small ( Steckcr 2007) . We have assumed the same form 
of the TT^-decay spectrum for starburst galaxies as for star-forming 
galaxies. 
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4.3. The Fermi Spectrum and Unresolved Sources vs. 
Truly Diffuse Mechanisms 

The Fermi observations have placed significant con- 
straints on extragalactic da.rk matter annihilation 
(iCireUi et al.ll2010l:lAb"do et al.|[20l0al: lAckermann et al.l 
l2010al) . Currently, there is no evidence of quark- 
annihilation features and spectral lines seen in the EGB 
spectrum, features that would be a clear annihilation sig- 
nal ( see e.g., Stec ker fc Tylka 1989a; Rudaz & Stecken 
119911) . The observed spectrum does not match that 
expected from dark matter annihilation, placing con- 
straints on any dark matter a nnihilation contribution to 
the EGB (lAbdo et al.ll2010aD . Therefore, it is probable 
that dark matter annihilation 7-rays, if present, provide 
only a minor contribution to the EGB. 

The same argument about matching spectra can 
be made regarding the contribution from electromag- 
netic cascades produced by very high and ultrahigh 
energy cosmic-ray interactions as the resulting spec- 
trum would be significantly harder tha n the observed 
spectrum (Kalashev et al. 2009: Bcrczinskv et al.l 120111 : 
lAhlers et alJI2010l: IVentersll2OlO0 . 

5. RESULTS 

The calculated spectrum of the unresolved FSRQ con- 
tribution to the EGB (see Sect ions 13 . 1 1 and I4?l1) is plotted 
in Figure [S] For c omparison, we incl ude the Fermi anal- 
ysis of the EGB (lAbdo et aLllMoil). two analyses^ of 
the EG RET EGB (ISreekumar et al.l 11991 IStrong et al.l 
l2004al) , and the calculation of the collective spectrum of 
unresolved FSRQs ignoring the effect of source confusion. 
Our results clearly show that the effect of source confu- 
sion is to reduce the number of resolved sources, increas- 
ing the collective intensity of unresolved blazars, particu- 
larly below ~ 1 GeV energy. Thus, accounting for source 
confusion modifies the predicted spectrum such that the 
EGRET and Fermi measurements of the EGB below ^ 1 
GeV are both compatible with unresolved FSRQs. In 
contrast, the better angular resolution of the Fermi-L AT 
above ^ 1 GeV allows it to resolve more blazars result- 
ing in a limiting flux that is dominated by the Fermi- 
LAT sensitivity rather than source confusion. Thus, the 
collective spectrum of FSRQs breaks at ^ 3 Ge\Q- At 
energies above ^ 1 GeV, the predicted collective spec- 
trum of FSRQs falls below the data points, though they 
are likely consistent with the data within the uncertain- 
ties in the galactic foreground emission model. Note also 
that the collective FS RQ spectrum exhib i ts mu ch less 
curvature than seen in IStecker fc SalamonI ^9M). This 
is because the spread in the S ID in our current model is 
much smaller than that of the IStecker fc SalamonI ()1996f ) 
model. 

In Figure m we plot the spectra of the unresolved star- 
forming galaxy contributions to the EGB calculated for 
the models discussed in Sections l3.2l and l42l For compar- 
ison, we include the spectrum of the unresolved starburst 
galaxy contribution alone that w e determined from th e 
best-fit IR luminosity function of IHopkins et al.l (j2010D . 

The two sets of EGRET data points result from two different 
estimations of the galactic foreground emission. 

The actual break should be more gradual since in our calcu- 
lations we used the approximate broken angular resolution curve 
shown in Figure |3] 




Fig. 5. — The collective spectrum of unresolved FSRQs. Solid 
green line: The spectrum accounting for source confusion. Dashed 
green line: The spectrum without accounting for source confu- 
sion. Black circles: The Fer mi mea surement of the spectrum 
of the EGB as determined in lAbdo et al.1 (2010i). Blue squares: 
The EGRET measurement of the spectrum of tfie EGB as deter- 
mined by Sreekumar et al. (1998) and confirmed by the analysis 
of Stecker et al. (2008) and S. D. Hunter (private communication). 
Red triangles: The EGRE T measur ement of the spectrum of the 
EGB as determined by Str ong et al.l l(2P04a). 

For the spectrum of starburst galaxies, we have assumed 
the same form of the tt*^ decay spectrum as for star- 
forming galaxies. The range in the calculations of the 
overall contribution to the EGB from unresolved star- 
forming galaxies spans about an order of magnitude in- 
dicating the degree of uncertainty in such a calculatiorF^ 
We note that even though our most extreme model 
could possibly explain the lowest energy Fermi data 
points (and possibly, within systematics, a couple oth- 
ers), it ca nnot explain the EG RET data points below 300 
MeV. The lStrong eTall (pOOla) EGRET data points (mi- 
nus the two highest energy data points) with the Fermi 
data points resemble a featureless power law, while the 
spectra of unresolved star-forming galaxies do not. No- 
tably, the data points show no indication of a Tr'^-decay 
"bump" at the energies at which the contribution of the 
star-forming galaxies should peak. 

6. DISCUSSION & CONCLUSIONS 

We have calculated the spectral shape of the contri- 
bution of unresolved FSRQs to the EGB assuming that 
the 7-ray luminosity of an FSRQ is, on average, pro- 
portional to its r adio luminosity (iGiroletti et al.l 120101: 
Abdo et al.ll2010il : iGhirlanda et al.ll2010l : iMahonv et al.l 
2010D . and also accounting for the effects of source con- 
fusion. We have demonstrated that the combination of 
the so urce density predicted by the Dunlop fc PeacocS 
(|1990l ) FSRQ radio luminosity function and the strong 
energy dependence of the Fermi-hAT angular resolution 
increases the contribution of unresolved FSRQs to the 
EGB at energies below 1 GeV. The resulting overall spec- 
trum predicted by the fit to the Fermi source count distri- 

Though, we note that each individual model is subject to its 
own uncertainty. As such, the degree of uncertainty is likely even 
more than an order of magnitude. Within the range of our various 
predictions of the EGB fro m star forming gala xies, we agree with 
the re sults of the model of [Fields et al.l ||201CI ) and IMakiva et all 
l(20Tll 'l. 
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Fig. 6. — The collective spectrum of unresolved star- forming 
galaxies. Dashed indigo line: Tlie spectr um deter mined from the 
strong coupling model (see Sections l3.2.3l and l4.2.3l Solid blue line: 
The spectru m dete rmi ned fro m the IR luminosity function model 
(see Sections 13.2.21 and I4.2.2")l . Dot-dashed yellow line: The spec- 
trum determined from the IR luminosity function model assuming 
no gas evolution. Dashed red line: The spectrum determined from 
the Schechter function model (see Sections l3.2.1l and l4.2.1ll . Double 
dot-dashed line: The spectrum of the starburst contribution alone 
determined from the IR luminosity function model. 

bution reproduces well the spectrum of the EGRET and 
Fermi EGB measurements below 1 GeV, but falls below 
the data points above 1 GeV. We have also calculated 
the spectral shape of the contribution of unresolved star- 
forming galaxies to the EGB for several relations for the 
7-ray luminosity of a star-forming galaxy. We find that, 
depending on the model, the overall amount of the contri- 
bution of star-forming galaxies to the EGB may be more 
or less significant, though regardless of the model con- 
sidered, the spectrum of unresolved star-forming galaxies 
is unable to explain the combined spectrum of the low- 
energy EGRET EGB measurements and the Fermi EGB 
measurements. Similar calculations for starburst galax- 
ies alone indicate that they account for at most about 1% 
of the EGB, in agreement with the conclusion reached by 
ISteckeiKpOOl . 

The similarity of the collective spectrum of unresolved 
FSRQs to the combined spectrum of the EGRET and 
Fermi EGB measurements, as demonstrated by our re- 
sults, is striking. In fact , we note that as predicted in 
IStecker fc SalamonI ()1999[ ) , the inclusion of the effect of 
source confusion in the calculation could provide an ex- 
planation for the similarity between the EGRET and 
Fermi EGB measurements at energies of hundreds of 
MeV. The density of FSRQs predicted by the model is 
sufficiently large such that at these energies, Fermi would 
not be able to resolve many more FSRQs than EGRET 
did, and the FSRQ contribution to the EGB would re- 
main the same for Fermi as for EGRET. Thus, if unre- 
solved FSRQs do comprise the bulk of the EGB emis- 
sion, then one would expect such similarity between the 
EGRET and Fermi measurements at these energies. At 
energies above 1 GeV, the Fermi-hAT angular resolution 
improves substantially with respect to that of EGRET, 
and as such, Fermi would be able to resolve more blazars 
at higher energies than EGRET could, resulting in a de- 
crease in the Fermi EGB with respect to the EGRET 



EGB, an effect which is possib ly indicated by comparing 
the EGRET and Fermi resultS 

In contrast, no such high-energy separation between 
the EGRET EGB and the Fermi EGB is predicted for 
star-forming galaxies as all but the closest are too faint 
to be resolvable by Fermi. Furthermore, we note that the 
EGRET EGB measurements provide no indication of a 
turnover in the spectrum as would be expected if un- 
resolved star-forming galaxies comprise the bulk of the 
EGB emission. Rather, the spectrum of unresolved star- 
forming galaxies is inconsistent with the combined spec- 
trum of the EGRET and Fermi EGB measurement J^. 
We also note that the lack of a turnover in the EGRET 
data is not simply the result of systematics (S. D. Hunter, 
private communication), since the uncertainties in all of 
the galactic foreground models used to determine the 
EGB from the EGRET and Fermi data are quite small 
at these energies. Finally, we note that at energies above 
~ 1 GeV, the spectrum of unresolved star-forming galax- 
ies is steeper than the spectra of the EGB data. As such, 
we conclude that however significant the contribution of 
star-forming galaxies to the EGB may be, it is not suffi- 
cient to explain the EGlfl. 

Within the range of our various predictions of the EGB 
from star forming gal axies, we agree wit h the results 
of the models of both [Fields et"aLl (|2010D (which sug- 
gests that star-forming galaxies may comprise the bulk 
of the EGB) and Makiya et al. (2011) (which suggests 
that star-forming galaxies can account for less than 10% 
of the EGB). This underscores the range of uncertainty 
in the calculation for star-forming galaxiej^. 

The featureless spectrum of the EGB deduced by 

A caveat is that the uncertainties in the subtraction of the 
galactic foreground emission at the higher energies are considerable 
owing to the uncertainty in the distributions of both gas and cosmic 
rays in the Galaxy. Furthermore, the instrumental backgrounds of 
EGRET and the Fermi-hAT are different, so it is difHcult to make a 
direct comparison between the two. We should also note that in our 
calculations, we have neglected the population of BL Lacs, which, 
due to their hard spectra, are likely to have more of a contribution 
at energies above ~ 10 GeV. Notably, Fermi has resolved as many 
BL Lacs as FSRQs. 

As previously noted, the Fermi-hAT was designed to reach 
its optimal effective area for 7-rays with energies near and above 
~ 1 GeV, whereas EGRET was designed to reach its optimal ef- 
fective area for 7-rays with energies near and above ^ 100 MeV. 
Also, the Fermi-hAT detector has a significantly higher instrumen- 
tal background at 100 MeV than EGRET did (S. D. Hunter, private 
communication). Thus, the EGB was not reported by Fermi for 
energies below 200 MeV IIAbdo et al.ll20l0il1 . 

The effect of Compton interactions mentioned in Section 13.21 
does not alter this conclusion as it only modifies the spectrum 
above 10 GeV for normal galaxies (Strong et al. 2010) and the 
starburst galaxy contribution to the EGB is negligible (See Figure 

El) 

One noteworthy difference between our model and that of 
[Fields et aLj (2010) is that in ord er to relate t he ga s mass of a 
galaxy to its star formation rate, ll<'ields et al.l (120101 ) makes use 
of the Schmidt-Kennicutt relation. However, in doing so, they 
estimate the disk sizes of galaxies to high redshifts. Given that 
the uncertainties in these quantities are likely to be considerable 
(and given the uncertainty already present in the calculation), we 
considered alternative approaches. Nevertheless, we note that the 
Schmidt-Kennicutt relation was included in the inputs to both the 
IR luminosity models and the Schechter function model. In the IR 
luminosity model, we tested the impact of changing the Schmidt- 
Kennicutt law by performing the calculation for the Hopkins et all 
(j2010) IR luminosity function calculated using a steeper Schmidt- 
Kennicutt relation and found that it had very little impact on our 
results. 



EGB Co 



Fermi is intriguing when one considers the possibil- 
ity of features that could arise from phenomena such 
as breaks in blazar spectra, absorption of high-energy 
7-rays from unresolved blazars, 7-ray emission from 
unresolved star-forming galaxies, 7-ray emission from 
dark matter annihilation, and 7-rays from electro- 
magnetic cascades initiated by very high and ultra- 
high energy particle interactions with the extragalac- 
tic background light. The spectra of these poten- 
tial contributions to the EGB differ considerabl y from 
that of the FSRQs (Silk & Srcdnicki 1984; Stccker etal] 
1985; Rudaz & Stccker 1988; S tecker fc T vlka 1989a^ 
Rudaz & Stcc ker. .1991,: ,Ullio et all 1 2002: Ando et al 



2007: IKalashev et al. 



200s ; iBerezinskv et al. 



20091: iSiegal-G askins & Pavlidou 
I20T1I Uhlers et al.ii2010l: iVenters 



2010f l. However, recent Fermi observations have 
placed significant constraints on dark matter annihilation 
ICircll i ct al., 2010.; ,Abdo et al.,, 2010a; Ackermann et al.l 
2010al) . and presently there is no clear evidence of anni- 
hilation features above the background continuum. As 
such, it appears that any putative contribution to the 
EGB from dark matter annihilation is relatively minor. 
The possible contribution to the EGB from electromag- 
netic cascades is constrained by the relative steepness of 
the EGB spectrum, though cascades co uld play a role at 



higher energies jKalashcv ct al. 200^ 



Berezinskv et al.l 



120111: lAhlers et aLri2010i: iVentersi |2dT3)- An apparent 



explanation for the featureless power-law spectrum of 
the EGB as presently deduced could be that unresolved 
blazars provide the dominant contribution to the EGB, 
given that their collective spectrum is roughly consistent 
with that of the EGB. 

Theref ore, we conclude t hat, contrary to the result 
given by lAbdo et al.l (|2010kl ) , the Fermi observations do 
not rule out the possibility that the EGB is dominated 
by emission from unresolved blazars. 
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